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The Kirkwood-Buff integrals for two-component mixtures in one-dimensional systems are cal- 
culated directly. The results are applied to square-well particles and found to agree with those 
obtained by the inversion of the Kirkwood-Buff theory of solutions. 



I. INTRODUCTION 

In Part I of this series^ one of us has studied both 
the global and the local properties of mixtures of sim- 
ple particles in one-dimensional system. This work has 
been part of a more general advocacy in favor of the 
study of local properties of liquid mixtures.— Instead of 
the traditional study of mixtures based on the global 
properties, such as excess Gibbs energy, entropy, volume, 
etc. we have advocated a shift in the paradigm towards 
focusing on the local properties of the same mixtures, 
such as affinities between two species (embodied in the 
Kirkwood-Buff integrals), and derived quantities such as 
local composition, preferential solvation, and solvation 
thermodynamic quantities. 

The local properties, though equivalent to and deriv- 
able from the global properties, offer a host of new in- 
formation on the local environments of each molecular 
species in the mixture. This information is not conspic- 
uous from the global properties. Therefore, the study of 
the local quantities offer a new and more detailed and 
interesting view of mixtures. 

In this paper we have recalculated the Kirkwood-Buff 
integrals (KBI) directly for two-component mixtures of 
particles interacting via square-well (SW) potential. 

In the next section, we outline the derivation of the 
pair correlation functions for two-component systems in 
ID system, for arbitrary nearest-neighbor interactions. 
In section UTIl we present a sample of results for mixtures 
of SW particles. It is shown that the results are in quan- 
titative agreement with those obtained in Part I, which 
were based on the partition function method and the in- 
version of the Kirkwood-Buff (KB) theory of solution. 
We have also calculated the limiting values of the KBI 
when one of the species has a vanishing mole fraction, 
which we could not have done from the partition func- 
tion methods. 

Another question examined both numerically and the- 
oretically is the deviations from symmetrical ideal solu- 
tions and its relation with the stability of the mixtures. 
It is shown that no miscibility gap can occur in such mix- 
tures. 



II. THEORETICAL BACKGROUND 

It is known that the correlation and thermodynamic 
properties of any one-dimensional homogeneous system 
in equilibrium can be derived exactly, provided that ev- 
ery particle interacts only with its nearest neighborship 
The aim of this section is to present a self-contained sum- 
mary of the exact solution. Although the scheme extends 
to any number of components,— here we focus on the two- 
component case. 

A. Correlation functions 

Let us consider a binary one-dimensional fluid mix- 
ture at temperature T, pressure P, and number densities 
Pa (ol = A,B). The particles are assumed to interact 
only between nearest neighbors via interaction potentials 
U a p(R). Before considering the pair correlation functions 
9ap{R)i it is convenient to introduce some probability 
distributions. 

Given a particle of species a at a certain position, let 
Pal(R)dR be the conditional probability of finding as its 
£th neighbor in some direction a particle of species (3 at 
a distance between R and R + dR. If I > 2 it is obvious 
that the (£ — l)th neighbor of a in the same direction 
(being located at some point R' between and R) is also 
a first neighbor of (J. Therefore, the following recurrence 
condition holds 

P%(R)= E [ R dR'p^(R')p^(R-R'), (2.1) 

■y=A,B J ° 

where P^p^R) is the nearest-neighbor probability dis- 
tribution function. On physical grounds^ the ratio 

Pal (R)/PaB C^) must become the same for a = A as 
for a = B in the limit of large R, i.e., 

ta -fffi . lim 4jiffl. (2 . 2) 

This relation will be used later on. The total probability 
density of finding a particle of species /3, given that a 



2 



particle of species a is at the origin, is 

oo 

p a p{R) = X)pg(iJ). (2.3) 

1=1 

The convolution structure of Eq. (|2.ip suggests the in- 
troduction of the Laplace transforms 

/>oo 

Jo 

P a p(s) = \ dRe- sR p a p(R), 
Jo 

so that Eq. (|2.ip becomes 



(2.4) 



7=A,B 



(2.5) 



Equation (|2.5[) allows us to express pTg(s) in terms of the 
nearest-neighbor distribution as 

pW( S )=[pW( S ) 



(2.6) 

where p^ (s) is the 2x2 matrix of elements pf^(s). F rom 
Eqs. ((231) and we get 

OO „ 

P( S ) = £ [p«( S ) 



(=i 



= p (1) (s)- fl-p (1) (s 



(2.7) 



where p(s) is the 2x2 matrix of elements p a f}{s) and I 
is the 2x2 unity matrix. 

Now, notice that the pair correlation function g a p(R) 
and the probability density p a p(R) are simply related by 
Pap{R) = Ppgaf3{R) or, equivalently in Laplace space, 



where 



/ dRe~ sR g af j(R) 
Jo 



(2. 



(2-9) 



is the Laplace transform of g a /}(R)- Therefore, thanks 
to the one-dimensional nature of the model and the re- 
striction to nearest-neighbor interactions, the knowledge 
of the nearest-neighbor distributions Vap{R) suffices to 
obtain the pair correlation functions g a p(R). More ex- 
plicitly, from Eqs. (|2.7|) and l|2.8|) the Laplace transforms 
9ap(s) are found to be 



9aa{s) 



9bb{s) 



1 

PT 



>AA 



(s) [l-Q BB (s)} + Q AB (s) 



xaD(s) 



(2-10) 



1 Qbb(s)[1-Qaa(s)]+Q 2 ab (s 



PT 



x B D{s) 



-, (2.11) 



9ab{s) 



Qab(s) 



p T s /x A ~x B ~D{sy 



(2.12) 



where pr = Pa + Pb is the total number density, x a = 
Pal Pt is the mole fraction of species a, and we have 
called 



*p( s ) 



(2.13) 



D(s) = [1 - Qaa(s)} [1 - Qbb(s)] - Q\ B {s). (2.14) 
The KBI in the one-dimensional case are defined by 



G a B — 2 



/•OO 

/ dR[g aP (R) - 1}. 
Jo 



(2.15) 



In terms of the Laplace transform g a f3(s), Eq. (|2.15p can 
be rewritten as 



G 



a/3 



2 lim 

s^O 



9a(3{s) 



(2.16) 



We see that only the nearest-neighbor distribution 
p^l(R) is needed to close the problem. It can be 
proven^, that p^p(R) is just proportional to the Boltz- 
mann factor e - u ^e>( R )/ k B T times a decaying exponen- 
tial e~^ R , where the damping coefficient is £ = P/k B T. 
Therefore, 



-U a!3 {R)/k B T e -tR 



(2.17) 



where the proportionality constants K a p = Kp a (which 
of course depend on the thermodynamic state of the mix- 
ture) will be determined below by applying physical con- 
sistency conditions. Taking Laplace transforms in Eq. 
(I2.17|) and inserting the result into Eq. (|2.13p we get 



lapis) = y/X a X/3K a 0Sl a p(s + £), 



where 



dRe- sR e- u ^' kBT 



(2.18) 



(2.19) 



is the Laplace transform of e~ Uafi ^ R ^ kBT . 

To recapitulate, given the interaction potentials 
U a p(R) and given a particular thermodynamic state 
(P,T, xa), the three correlation functions are obtained 
(in Laplace space) from Eqs. (12.10j) - (|2.12[) . supplemented 
by Eqs. (j2~Tl) . (|235|l . and (l2~Tgi) . 



B. Equation of state 

In order to close the exact solution, it only remains to 
determine the total density p? (equation of state) and the 
amplitudes K a p as functions of P, T, and x A = 1 — x B . 
As said above, they can be easily obtained by applying 
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basic physical conditions. First, note that Eq. (|2.2|) es- 
tablishes the following relationship 



K 2 AB = K AA K 



BB- 



(2.20) 



Next, the physical condition lim^-joo g a p(R) = 1 implies 
that g a p(s) — > 1/s for small s. According to Eqs. (|2. 10|) 
(|2~T2|) . this is only possible if D(Q) = 0, so that D(s) -> 
D'(0)s for small s, where D'(s) — dD(s)/ds. Thus, one 
has 

[1 - Qaa(0)}[1 - Qbb(0)} - Q 2 AB (0) = 0, (2.21) 



Pt 



Qab(0) 



Pt 



Pt 



Qaa{0)[1-Qbb( 



QabW 



x A D'(0) 



Qbb(0)[1 - Qaa{0)} + Q 2 AB (0) 
x B D'(0) 



(2.22) 



(2.23) 



(2.24) 



Elimination of pt between Eqs. (|2.22p ~ (|2.24[) yields two 
coupled equations which, together with Eq. (|2.2ip . gives 



K AA 



1 - X B K AB fl AB (£) 

x A tt AA {£,) 



(2.25) 



A" 



BB 



i - x A K AB n AB ((,) 

x B ^t BB {Cj 



(2.26) 



Insertion of Eqs. fUg) and (|2~2"6"1) into Eq. ([2~2TJj) allows 
one to obtain a quadratic equation for K AB whose phys- 
ical root is 



K 1 1 — \/l — 4x A x B (l — R) 

AB ~M() 2^ B (l-A) ' (2 ' 27) 



where we have called 



R 



= n AA (fln BB (fl 



(2.28) 



It is interesting to note that, since K a p and f2 Q ^ 
are positive definite, Eq. (|2.25[) and (|2.26[) imply that 
x a K AB Vt AB {£) < 1 for a = A,B, i.e., 



K AB n AB (0 < min (—,— )< 2. (2.29) 

\x A x B 



Finally, the density pt is obtained from either of Eqs. 
([2~^ - (pT24]) . The result is 



Pt(P,T, x a ) 



x A K AA n' AA (t) + x 2 b k bb q' bb {o + 2 XA x B K AB n AB (t) ' 



(2.30) 



where Q' a n(s) is the first derivative of £l a p(s). 

Equations (pT^ - IpE?^ and (pOOj) complete the full 
determination of g a /3(s) and the equation of state for 
any choice of the nearest-neighbor interaction potentials 
U a p{x) and of the thermodynamic state (P,T,x A ). 



C. Kirkwood— Buff integrals 



The KBI G a /3 can be derived, according to Eq. (|2.16p . 
by expanding sg a j3{s) in powers of s as sg a p(s) = 1 + 
\G a ps + ■ ■ ■ and identifying the linear term. After some 
algebra one gets 



Gab = PtJ + 2 



"WO 
Om(0' 



(2.31) 



Gbb — PtJ — 2 



x A K AA n' AA (i) 2 



x B K AB tt AB (£) p T x B ' 
where 

+2a; A a; B A Ai 30 AB (e) 



(2.33) 



(2.34) 



The knowledge of the KBI allows us to obtain the (re- 
duced) isothermal compressibility 



X = k B T 



dpT 
dP 



(2.35) 



T,x A 



by means of 



G AA = PtJ - 2 



x B K BB n' BB (Q 2_ 

x A K AB Q, AB (£) p T x A ' 



(2.32) 



1 + p T x A x B A AB 
+p T x A x B (G aa G bb — G AB j 



[1 + p T {x a G aa + x B G BB ) 

2 ^ (2.36) 
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where 



-0.5 



A AB = Gaa + Gbb — 2Gai 



(2.37) 



It can be checked that the resulting expression of \ 
(which, due to its length, will be omitted here) coincides 
with the one obtained as x = {®Pt/&£)tx f rom Eq. 
(I2.30p . This confirms the exact character of the solution. 

Making use of Eqs. (|2.30[) - (|2.33p . it is easy to prove 
that 

1 + ptx a x b A ab = v J. !e , -l, (2.38) 
Kab^ab(U 

which, according to Eq. (|2.29p . is a positive definite quan- 
tity. More explicitly, from Eq. (|2.27[) we have 

1 + PTXAXBAAB = 4x A X B (l - R). (2.39) 

Therefore, the denominator in Eq. (|2.36|) never vanishes 
and the isothermal compressibility is well defined. This 
agrees with van Hove's classical proofi that no phase 
transition can exist in this class of nearest-neighbor one- 
dimensional models. 

Let us now obtain the KBI in the infinite dilution limit 
x A -> 0. In that limit, Eqs. |2~25)) - (i2~27|) and (f2~30"l) 
become 

1 1 

K A b = 



K AA = 775 — 777, i^BB 



&ab(0' 



Pt 



n BB (Q 

V'bb(0' 



(2.40) 
(2.41) 



Analogously, from Eqs. (|2.31[) - (|2.34p one gets 

c Aa = "HQ I 2 n ' ABi0 



Gaa 



^(0 , A n' AB (0 



G 



BB 



A AB = 2n' BB (0 



<W0 
i 



-2 



& AB (0 



'^bb{0' 



(2.42) 
(2.43) 
(2.44) 
(2.45) 



n BB (0 n AB (0. 

Note that special care is needed to obtain K AA and G AA . 

D. Chemical potentials and solvation Gibbs 
energies 

Finally, let us get an explicit expression for the chemi- 
cal potential. From the KB theory of solution we have^ 

1 / dp A \ 1 p T x A A AB 

k B T \dx A J PT x A 

1 

x A 



1 + p T x A x B A AB 




FIG. 1: The KBI Gap for hard rods of different diameters 
<Jbb/<Jaa — 2 and Pa A A/kBT = 1. 



where in the last step we have made use of Eq. (|2.39[) . 
Integration over x A yields 



PA 
k B T 



const + \nx A + In [1 — 2x B (l — R) 



+ ^l-Ax A x B (l- R) 
For pure A (x B — 0), we have 

const + In 2. 



Va 
kuT 



(2.47) 



(2.48) 



The solvation Gibbs energy of A in pure A may be ob- 
tained from (|2.47j) as£& 



Ap* A = p A - k B T\n{p A A A ), 



(2.49) 



where A A = h/-\/2irm A k B T is the momentum partition 
function of A in one-dimensional systems. Similarly, 



a Ma P = Ma - k B Tln(p A A A ), 



(2.50) 



where p p A is the density of pure A at the same T and P 
as the mixture. Taking the limit x B — > in Eqs. (|2.25p 
and (|2.30[) one has 



P P A 



(2.51) 



The excess solvation Gibbs energy relative to the sol- 
vation Gibbs energy in pure A is defined as 

AAp* A = Ap* A - Ap* A (2.52) 



This quantity may be calculated from (|2.47|) - ()2.52|) with 
the result 



1 y/l-Ax A x B {\-R)-l k B T 
x B ^1- Ax A x B {\- R) 

(2.46) 



In 



1 



--x B {l-R) + -yJ\-Ax A x B {\-R) 



ln^. 

Pt 



(2.53) 
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FIG. 2: The KBI Gaa for SW particles with parameters given 
in and k B T/\e AA \ = 1, Pa/k B T = 1. The line s are 

obtained from the exact expressions presented in Sec. Ill CI 
while the circles are the data obtained in Ref. Q]. 
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0.0 0.2 0.4 0.6 0.8 1.0 

X A 

FIG. 3: The KBI G BB for SW particles with parameters given 
in (J3T7J| and k B T/\e AA \ = 1, Pa/k B T = 1. The line s are 
obtained from the exact expressions presented in Sec. Ill CI 
while the circles are the data obtained in Ref. Q]. 



III. A SAMPLE OF RESULTS 

Let us start considering a binary system composed of 
(additive) hard rods of different diameters (lengths) <jaa, 
o be, and gab = (&aa + cbb)/2- The Laplace function 
fl a0 (s) denned by Eq. (|2H?j) is 

fia/s(a) = ■ (3-1) 

s 

In this case the parameter defined in Eq. (|2.28[) is R = 1 
and thus the limit R — ► 1 must be taken in Eq. (|2.27[) 
with the result Kab — The general scheme 

of section [IT] can be used to obtain the KBI explicitly: 

„ OAA + OBB +&AAOBB lo ^ 

Gab = ~ , ; r, (3.2) 

1 + i{XA(JAA + XB&BB ) 



Gbb = Gab + oaa - &bb, (3.4) 

so that Aab = 0. Figure Q] shows the values of G a p for a 
diameter ratio <jbbI&aa = 2 and a thermodynamic state 
PcTAA/ksT = 1. These results are in perfect agreement 
with those calculated in Part ii 

Having established that the programs give the correct 
results for hard rods, we next present results for a mix- 
ture of particles' interaction via SW potential of the form 

{OO, R<(T a p, 
e Q /3, cr Q/3 < R < <r a/3 + 6 a p, (3.5) 
0, R < a a p + 5 aj3 . 



Gaa = Gab + &bb ~ vaa. 



(3.3) 



where e a p < 0. For this SW potential the Laplace func- 
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0.0 0.2 0.4 0.6 0.8 1.0 

X A 

FIG. 4: The KBI Gab for SW particles with parameters given 
in (J3T7J| and k B T/\e A A\ = 1, Pa/k B T = 1. The line s are 
obtained from the exact expressions presented in Sec. Ill CI 
while the circles are the data obtained in Ref. [l|. 



s=-l 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 5: Values of Aab for SW particles with parameters 
given in and k B T/\e A A\ = 1, Po/k B T = 1. The lines 

are obtained from the exact expressions presented in Sec. Ill Cl 
while the circles are the data obtained in Ref. Q]. 



tion £l a p(s) is 



f2 Q/ 3(s) 



-e a /3/k B T 



(^ e -e ag /k B T _ ^ 



,sS n 



(3.6) 

and again the general results of section [U provide the 
KBI explicitly. 

We have taken the following values for the potential 
parameters: 



OAA = f BB = OAB = CT, 

Saa = Sbb = Sab = -<r, 
5 

ess 



(3.7) 



SAA 



£, £AS — — \JtAAtBB- 



The thermodynamic variables are T, P, and xa- In 
all the calculations we choose ksT I\^aa\ = 1 and 



Pa/kBT = 1 to compare the present results with those 
of Part I. 

Figures OB] show the values of Gaa, Gbb, and Gab 
for these systems for various values of e ranging from 
e = -0.001 to e = -1, and from e = -1 to e = -113 
Figure [5] shows the values of Aab = Gaa + Gbb— 2Gab 
in the entire range of composition. In all the cases the 
agreement with the results of Part I is quantitative. 

The KBI in the infinite dilution limit (xa — > 0), as 
obtained from Eqs. (j2l5|) - (|2l4|) . are plotted in Fig. [6] 
as functions of — e for the same system as that of Figs. 
[2HS1 We observe that both Gab and Gbb arc hardly 
sensitive to the value of e. In contrast, the solute-solute 
KBI, Gaa, is strongly influenced by the solvent-solvent 
potential depth, increasing both for small and for large 
values of | e | . A careful inspection of the explicit expres- 
sions (|2.42p - (|2.44p in the limit |e| — > oo shows that, while 
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FIG. 6: The KBI Gap in the infinite dilution limit (xa — * 
0) for SW particles with parameters given in (|3.7p and 
k B T/\e AA \ = 1, Pa/k B T = l. 



Gab and Gbb tend to the same constant value, G aa di- 
verges as G A a ~ exp [(|e BB | - 2|cab|) /fc B T]. This phe- 
nomenon might be relevant to the study of hydrophobic 
interactions, as discussed in Ref. U 



IV. DISCUSSION AND CONCLUSION 

In Part I we calculated all the KBI in an indirect wayi 
We first calculated the excess functions from the partition 
function of the system, then we used the inversion of the 
KB theory^ to calculate the KBI. This lengthy procedure 
might have introduced accumulated errors. Some readers 
of Part I have expressed doubts regarding the reliability 
of the results calculated along this procedure. In fact 
some have also claimed that there might be a miscibility 



gap which we might have missed by this indirect and 
lengthy calculations. 

In this paper we have repeated the calculations of the 
KBI directly, from the same program that was designed 
to calculate the pair correlation functions in mixtures of 
two components in ID system. 

The agreement between the two methods was satisfy- 
ing, it also lent credibility to the inversion procedure and 
encouraged us to extend the calculations of the KBI for 
aqueous like mixtures £ We hope to report on that in the 
near future. 

Regarding the question of miscibility gap we have 
shown that the inequality 

1 + ptXaxb^ab > (4.1) 

always holds in these mixtures, as shown by Eq. (|2.39D . 
From the KB theory^ of solution we have the equation 

(d 2 g\ = J_ /dfiA\ = ksT 

\dx 2 A ) PT XB\dxAj PT xaxb(1 + ptxaxbAab)' 

(4.2) 

where g = G/(Na + Nb) is the Gibbs energy of the 
system per mole of mixture. It follows from (|4.1[) and 
(|4.2p that g is everywhere a concave (downward) function 
of xa- Therefore, there exists no region of compositions 
where the system is not stable, hence no phase transition 
in such a system. 
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